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Abstract: We describe a simple analytical method for effective summation of series, 

including divergent series. The method is based on self-similar approximation theory 
resulting in self-similar root approximants. The method is shown to be general and applicable 
to different problems, as is illustrated by a number of examples. The accuracy of the method 
is not worse, and in many cases better, than that of Fade approximants, when the latter can 
be defined. 

Keywords: Asymptotic series; Effective summation; Root approximants; Self-similar 
approximation theory 


1. Introduction 

In numerous cases of applied mathematics and mathematical physics the solutions to problems can 
only be represented as series derived by means of some kind of perturbation theory or iterative procedure. 
A great majority of such series is even divergent, having meaning only as asymptotic series for an 
infinitesimally small expansion variable. While the considered problems often require to consider finite 
values of this variable, sometimes even very large values. The standard way of treating such asymptotic 
series, for the purpose of their extrapolation to the finite values of the variable, is by invoking the Fade 
approximants [1]. The latter, however, exhibit several deficiencies limiting their applicability, as is 
discussed in Refs. [1,2], for instance, such a notorious deficiency as the appearance of spurious poles. 
Another weak point is the ambiguity of choosing one of the Fade approximants Pm/n from the table of 
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many admissible, for each series of order k, variants satisfying the condition M + N + 1 = k. Also, in 
the limit of a large variable x the approximant Pm/N ix) behaves as x^~^. Hence, only integer powers 
of X are allowed. It is possible to improve the results by employing the modified Fade approximants [3], 
corresponding to the power with choosing the appropriate value of 7 satisfying the large-variable 

limit. 

In the present paper, we show that it is possible to formulate a general method for effectively 
extrapolating and interpolating asymptotic series. The method enjoys the following advantages: (i) It is 
unambiguously defined for each given series of order fc; (ii) It allows for the treatment of large-variable 
behavior of any type, whether with integer, rational, or irrational powers; (iii) Being more general, it is 
not less accurate than the method of the Fade approximants, when the latter exist, in many cases, being 
more accurate. 

In the great majority of realistic situations, only a few terms of asymptotic expansions are available. 
Therefore, in the examples below, we do not consider very large series, showing that even several of 
terms allow us to derive quite accurate approximations. 


2. Self-similar root approximants 


Suppose we are interested in finding a real function fix) of a real variable x. However, this function 
is defined by a complicated equation that cannot be solved exactly. But, applying a kind of perturbation 
theory, we can derive the small-variable behavior of this function 


fix) ^ fkix) (x ^ 0 ) , 

represented by asymptotic series, with the fc-th order expansion 


fkix) = foix) ( 1 + 


a„x 


n=l 


where 


foix) = . 

Sometimes, the large-variable behavior of the function 

fix) ~ f^^\x) ix —)■ 00 ) 

is also known and can be represented by an expansion over 1 /x as 

f^^\x) = f^ix) + ^ 


n=l 


with 


/oo(x) = Bx^ . 


( 1 ) 


( 2 ) 


( 3 ) 


( 4 ) 


( 5 ) 


( 6 ) 


For what follows, it is convenient to deal with the ratio fix)/foix), which at small variable x ^ 0 
behaves as 

fix) fkix) 


foix) foix) 


=1+5: 


a^x 


n=l 


(7) 
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and at large values of the variable x —oo it tends to 

/o(x) /o(x) A 


( 8 ) 


The extrapolation of the small-variable expansions to the large-variable region ean be done by 
means of self-similar approximation theory [4-8]. In this approaeh, the transfer from a fc-th order 
approximation, say, a small-variable expansion, to the higher orders of approximations are treated as 
the motion with respeet to the approximation order k playing the role of diserete time. Construeting 
a dynamieal system, whose trajeetory is bijeetive to the sequenee of approximations, makes it feasible 
to find a fixed point representing the sought funetion. The eonvergenee to the fixed point is governed 
by eontrol funetions. The self-similar approximation theory eombines the methods of optimal eontrol 
theory, dynamieal theory, and renormalization-group approaeh. We shall not go into the details and 
mathematieal justifieation of the self-similar approximation theory that has been thoroughly expounded 
in Refs. [4-8], but we shall use some of its eonsequenees. 

Employing this theory for the purpose of interpolation between the small-variable and large-variable 
regions, it is possible to eome [9-11] to the self-similar root approximant 


^ =((... (1 + A,xr^ + A2xY" + • • • + • (9) 

A theorem has been proved [12] stating that all parameters Ai and powers rii of approximant (9) are 
uniquely defined through the large-variable form (5). 

However, the root approximant (9) eannot be uniquely defined through the small-variable expansion 
(2). This hinders the applieability of the approximant (9), sinee in the majority of oases, the 
small-variable expansion is better known, providing a number of terms, while the knowledge of the 
large-variable behavior is limited by just a single term ( 6 ), often even without preeise data for the 
amplitude B. In order to extend the applieability of approximant (9) to be uniquely defined through 
the small-variable expansion, it is neoessary to impose some oonstraints on the powers rij. Suoh a 
straightforward oonstraint is the requirement that all parameters Aj of approximant (9) be involved in 
the definition of the large-variable limit, whieh implies the relation 

= ^ (j = l,2 ,...,fc-l), (10) 

with Uk = P — a. By expanding Eq. (9) in powers of x, it is easy to prove that all parameters Aj are 
uniquely defined through the eoeffieients aj of small-variable expansion (2). In addition, we ean require 
the validity of the limiting form ( 6 ), whieh improves aeeuraey. 

The self-similar root approximant (9), with eonditions ( 6 ) and (10), whose parameters Aj are uniquely 
defined by the aeeuraey-through-order proeedure and are expressed through the eoeffieients aj of the 
small-variable expansion (2), ean be ealled, for short, the root approximant. In the following seetions, 
we demonstrate that this root approximant provides quite aeeurate approximations for different problems, 
uniformly extrapolating the small-variable expansion (2), valid for x ^ 0, to the whole region of x G 
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3. Illustration by simple examples 

Before going to more eomplieated problems, we show the effieieney of the method by simple eases. 


3.1. Hard-core scattering problem 


Let us start the illustration of the method from the problem eonsidered by Baker and Gammel [3]. 
When ealeulating the seattering length of a repulsive square-well potential, one meets the integral 






dt , 


whose limit, as a; ^ cxo, equals vr/lS. Baker and Gammel state that this integral eannot be eorreetly 
evaluated by the standard Fade method. To solve the problem, they suggest a modified method employing 
a power of the Fade approximant. We show below that sueh integrals ean easily be treated by means of 
the root approximants. 

The small-variable expansion of this integral reads as 

/y»13 

- g - - 297675 ^ 5893965 ~ 166080925 ^ 10672286625 ‘ 

Comparing this with form (2), we have S'o(a:) = x/9. Sinee expansion (2) is in powers of we 
eonstruet the root approximants (9) using as a variable. Thus, the root approximant of third order is 

where the parameters are 


Til = 0.133333 , 712 = 0.012952, 713 = 0.016907. 

To fourth order, 

Slix) 

where 

Til = 0.133333 , A 2 = 0.012952 , 7 I 3 = 0.002757 , 7 I 4 = 0.004636 . 

To fifth order, 

^ + ^ 33 ;®^ +A 4 X^^ AA^x^^ 

where 

Til = 0.133333 , 7 I 2 = 0.012952 , 7 I 3 = 0.002757 , 



^ ((^(^{1 + Aix'^y + A2 x‘^'^ ^ + Asx^^ +A4X 


-1/8 


= 0.000578 , 


7 I 5 = 0.001285 . 
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And to sixth order, 

^ ^ ^ +^ 432 ®^ + 

where 

Ai = 0.133333 , A 2 = 0.012952 , A 3 = 0.002757 , 

A 4 = 0.000578 , As = 0.000137 , Ag = 0.000356 . 

All these approximants eonverge to 7r/15, as a: ^ cxo. The higher the approximant order, the faster the 
eonvergenee. 


3.2. Debye function 


The n-th order Debye funetion is defined [13] through the integral representation 

rx 


D{n, x) = 


n 


x" 


'0 


e* — 1 


dt . 


For |a;| < 27r and n > 1, it possesses the expansion 


D{n, x) ~ 1 


n 


X 


+ 




2fc ^2k 


2 (n + l) ^ ( 2 A; + n)( 2 fc)! 

in whieh B 2 k are Bernoulli numbers. At large x and Re n > 0, one has 

a. 


D{n,x) ~ 


x" 


(x —)■ cx), Re n > 0) , 


where 


Cn = nT{n + 1)C(?^ + 1) • 

Below, we eonsider the ease of n = 3, eorresponding to the Debye funetion 

D{x) = D{3,x) = — I — -r dt . 


x-^ 


e* — 1 


The small-variable expansion for the latter takes the form 


O ^ 

D{x) ~ 1 - o ^ a 2 kx‘^’" (x -)■ 0) , 


k=l 


in whieh 


B 


0 ‘ 2 k — 


2 k 


{ 2 k + 3){2k)\ 

While the large-variable behavior is given by the expression 


D{x) ~ ^ 


X"' 


(x —)> 00 ) , 
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with 



5 


19.481818 . 


Constructing the root approximant 



■ 3/5 


X 


9 


we eompare it with the exaet numerieal values of the function D{x) and find that D^{x) approximates 
well this function in the whole region of x G [0, oo], with the maximal error of 15% at x = 5. The best 
two-point Fade approximant of the same order, Pi/ 4 (x), is less aeeurate, yielding the maximal error of 
33% at X = 15. 

3.3. Fermi-Dirac integral 

The general form of the j-th order Fermi-Dirac integral is 



Its asymptotic expansions are known [14]. 

For concreteness, let us consider the zero-order ease that reduees to the function 


F(x) = F( 0 , x) = In (1 -f e^) . 


At small X, this function tends to In 2, and at large x, we have 


F(x) ~ X (x oo) . 


The root approximant 


-f Aixf + Asx^) 



1/5 


(x) = In 2 


3/2 


where 


Ai = 0.721348 , A 2 = 0.360674 , A 3 = 0.390257 , 
A 4 = 0.410334 , A 5 = 4.294519 , 


provides an accurate approximation for the funetion F (x) in the whole region of x E [0, 00 ] , the maximal 
error being 5%. The two-point Fade approximant P 3 / 2 (x) is slightly less aeeurate, with the maximal error 
of 6 %. 
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3.4. Fekete-Szego problem 

The problem of maximizing the absolute value of a funetional in subclasses of normalized functions 
is called the Fekete-Szego problem [15,16]. The Fekete-Szego functional is bounded by the function 



where 0 < x < 1. 

In order to consider the interval [0, cxo], as in other examples, we can use the change of the variable 



X 


1 — X 


Then z —)> oo as x —1. Expanding F{z) at small z gives 


Fiz) = f(x(z)) ^3-4z + 4z^- -z^ + - z^ - — z^ . 

K J JK K jj 3 3 


The root approximant F^{z) uniformly approximates the function F{z) on the interval z G [0, oo], with 
the maximal error about 10%. The two-point Fade approximant P 2 / 2 {z) is worse, having the maximal 
error twice larger than the root approximant F^{z). 

4. Some useful tricks 

It is important to mention some tricks allowing for the convenient use of the method. Below, we 
discuss the interchange of small-variable and large -variable limits and the problem of dealing with 
logarithms. 

4.L Inversion of expansions 

In the above examples, we have considered functions, whose expansions are better known for the 
small-variable limit, while a few, or just a single term, are available in the large-variable limit. But 
generally, the small-variable and large-variable limits are interchangeable. In those cases, when the 
large-variable expansion in powers oil/x provides a number of terms and this expansion enjoys better 
convergence properties, it is possible to inverse the small-variable limit to the large-variable limit by 
using the variable change x = 1/t. Then, instead of the function f{x), we consider the function 



t = — 

X 


1 


( 11 ) 


The small-variable limit (1) becomes the large-variable limit 



( 12 ) 


in which 



( 13 ) 
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with 


Fooit) = fo[-]=At- 


Conversely, the large-variable behavior (4) transforms to the small-variable behavior 


in whieh 


where 


F(t) ~ F,(t) = /M Q (( ^ 0) , 


Fp(i) =-FoW + 

n=l 


(14) 


(15) 


Fo{t) = U 


= Bt-f^ . 


After this ehange of the variable it is straightforward to employ the same proeedure of eonstrueting 
the root approximants, as is explained in See. 2. 

More generally, it is possible to use the ehange of the variable f = 1/x®, with a positive power s > 0, 
so that again t ^ oo, when x —)■ 0. 

4.2. Example of inversion 

As an illustration of the inversion proeedure, we give below a typieal example, diseussing it rather 
briefly, sinee the whole method of eonstrueting the root approximants is the same as before. 

Let us eonsider the partition funetion of the so-ealled zero-dimensional oseillator, or the generating 
funetional of zero-dimensional field theory, whieh is defined through the integral 

1 f°° 

I{x) = —;= / exp (— dip , 

TT / — 


where x plays the role of a eoupling parameter. In the weak-eoupling limit, one has [17] the asymptotie 
expansion 

OO 

/(x) ~ 1 -f anX" (x 0) , 

n=\ 

(-1)'" f 1 

On = r 2n + - 
/tt n! \ 2 


in whieh the eoeffieients are 


For instanee 


105 


ai = 


4 ’ 


0,2 — 


32 ’ 


0-3 — 


3465 
128 ’ 


and so on. 

The strong-eoupling expansion reads as 


/(x) ~ 1.022765 x-F^ - 0.345684 x'^/^ + 0.127846 x"®/^ {x ^ oo) . 


Here the strong-eoupling expansion provides a number of terms. Moreover, the absolute values of 
the eoeffieients in this expansion diminish with inereasing order, eontrary to the eoeffieients a„ in the 
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weak-coupling expansion, which grow as rf with increasing order n. This makes the strong-coupling 
expansion more suitable for constructing root approximants. 

Resorting to the change of the variable x = we consider the function J{t) = and follow 

the scheme of the previous section. We define the root approximants Jl{t) that give us the approximants 
/^(x) = for the sought function. Found in that way approximant has the maximal 

error of 5% for the whole range of x G [0, cxo]. For comparison, the Fade approximant Pi/ 2 ix) has the 
maximal error of about 20%, which is much less accurate. 


4.3. Dealing with logarithms 


It is worth paying attention to the problem of series involving logarithms, which often appear in 
physics applications. Such series do not yield any complication for the method of root approximants 
described here.There are two equivalent ways of treating such series. Thus, if a series contains the terms 
with x”, and with x"^ Inx, then it is admissible to consider as the terms of one order either those 
containing x"' and x"' In x or the terms and x"' In x. 

As an illustration, let us consider, e.g., the typical form of such a series involving logarithms as that 
one arising in the Nambu-Iona Lasinio model [18] and leading to the function 


/(x) = X 


Vi + 




X 


In 


1 + \/l + x^ 


X 


where x plays the role of mass. At asymptotically small x, it follows 


/(x) ~ X + 


— In 2 + In X X 


(x ^ 0) . 


While at large x, one has 




28x4 


X —)■ CXO 


Keeping in mind the dependence of the last expansion on 1/x^, it is convenient to use the variable 
2 ; = 1/x^. The root approximant, satisfying the required limits, has the form 


/I( 2 ^) — 2 ^((1 + ^ 1 ^) +^ 2 ^^) + \n{l-\-z) + , 

with all parameters uniquely defined by the given expansions. This expression approximates well the 
initial function /(x), with the maximal error of 2% at x ~ 2. Contrary to this, the best Fade approximant 
of the same order has the error of 11% at x ~ 1.5. 


5. Ground-state energy of electron gas 

Important and not trivial problems arise when studying the properties of charged systems [19-21]. 
Here we show how our method works for the case of homogeneous electron systems. 
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5.1. One-dimensional electron gas 

The Hartree-Fock part of the uniform electron energy is well known. The problem arises in calculating 
the correlation energy. The latter is usually presented in a reduced dimensionless form e{rs) as a function 
of the Seitz radius r^. High-density expansion for one-dimensional uniform electron gas [22] corresponds 
to small Vs, when for the correlation energy one has 

e{rs) ~ C -I- 0.00845rs (r^ —)■ 0) , 

where 

TT^ 

C = -= -0.027416 . 

360 

The low-density expansion [22] implies large when 

, . bi b2 . . 

sirs) ^ ^ (r. ^ cx)) , 

s Ts 

where 

bi = - “ i) ^ -0.168939 , 62 = 0.359933 . 

The root approximant, enjoying the same expansions, but valid for arbitrary reads as 

with the parameters 

Til = 0.493150 , TI 2 = 0.056122 , TI 3 = 0.004274 . 

Comparing the prediction of the root approximant with the data from diffusion Monte Carlo calculations 
[22] in the interval 0 < < 20, we find that the maximal error of £3 is 8 %. Fade approximants 

give the errors between 2% and 10%. Thus, Pi/ 2 (\A^) has the error of 2%, while To/3(\A7) has the 
maximal error of 10%. The Cioslowski interpolation method [23] results [22] in a better accuracy of 
1%. However, this method includes an additional parameter that is fitted from numerical Monte Carlo 
calculations. While our aim has been in constructing good approximations without fitting parameters, 
being based only on asymptotic expansions. The principal importance of avoiding fitting parameters is 
crucial for those problems where no exact numerical data are available. 


5.2. Two-dimensional electron gas 

Correlation energy of a homogenous two-dimensional electron gas was studied in several articles, 
e.g., in Refs. [24-31]. In high-density limit (small r^), the ground-state energy reads [31] as 

Eoivs) ~ ^ -f — -f e(rs) (r^ 0) , 

rj Ts 

where the first two terms constitute the Hartree-Fock energy, with 



4x/2 

dvr 


c-i = - 
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And the last term is the eorrelation energy 

- Co + c\rs Inr^ (r^ 0) , 

with the eoeffieients 

Co = -0.192495 , c[ = -s/2 = -0-0863136 . 

In the low-density limit (large r^) the asymptotie expansion for the eorrelation energy ean be written 
[26] as 

/ X I ^2 ^3 / ^ X 

c r, ^ — + ^ + — r, ^ cx) , 
rs r7 

where 

= -0.472189 , 62 = 0.4964 , 63 = 0.5297 . 

For intermediate r^, there have been suggested [27,28,30] several phenomenologieal expressions with 
parameters fitted from Monte Carlo ealeulations. Thus, Gori-Giorgi et al. [28] suggested the form 


^(r^) = Ao + {B^rs + Cor] + Dorfj In ( 1 + ^ ^ 3/2 


Eor, + + G^rl + 


with the parameters 


Ao = -0.1925 , Bi = 0.0863136 , Go = 0.057234 , Dq = 0.003362896 . 

Eq = 1.0022 , Fq = -0.02069 , Go = 0.34 , Hq = 0.01747 . 

This expression ean be used as a numerieal result for estimating the aeeuraey of approximate analytie 
formulas. 

The root approximant, satisfying all asymptotie expansions reads as 


= 


hi 


o \ 3/2 

Al \ A2 \ A3 / 1 

l + ^ +— + — ln l + 


r.. 


A4 A5 

j- 1 . j-2. 

' 3/2 ' J.2 


-1/2 


where the parameters are 


= -0.472189 , Al = 0.700849 , A 2 = 2.723702 , 

A 3 = 10.792193 , A 4 = -5.764339 , A 5 = 6.017150 , 
The error of this approximant is about 5%. 


6. Systems with spherical symmetry 

Finite quantum systems often enjoy spherieal symmetry. Below, we eonsider two examples of sueh 
systems that are important for applieations. 
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6.1. Energy of harmonium atoms 

An A^-electron harmonium atom is described by the Hamiltonian 

TV TV 

H = -y" (-+ a;Vf) + - V — 

2 ^ ^ ^ 2 ^ Vi, 

1 = 1 lytj 

where dimensionless units are employed and 


Vi = Fi 


nj = |ri - r^l 


This Hamiltonian provides a rather realistic modeling of trapped ions, quantum dots, and some other 
finite systems, such as atomic nuclei and metallic grains [32]. This is why the energy of harmonium 
atoms has been intensively studied [33-37]. Here we show that root approximants give a good 
approximation for the energy of such systems. We consider the ground-state energy of a two-electron 
harmonium. 

At a shallow harmonic potential, the energy can be expanded [23] in powers of to, so that 

E{u}) ~ Ekioj) (cu ^ 0) , 


with the truncated series 


For instance, to third order, we get 




Cr,U 


(2-rn)/3 


n=0 


E^[uj) — -f CiUJ -f C2UJ 


4/3 


with the coefficients 
3 


Co = 


24/3 


= 1.19055 , 


ci = - (3 + ^3) = 2.36603 , 


And for a rigid potential, the energy is approximated [23] as 

E{oi}) ~ E^^\uj) {oj 00 ) 


where 


To fourth order, one has 


where 


E^P\0J) = 


u 


(2-n)/2 


n=0 


7 

C2 = — 2-2/3 = 0.122492 . 
36 


= bou + + ^2 + bsoj ^ 


60 = 3 , 61 = W-= 0.797885 , 62 = -- fl--+ln2 )=-0.077891 , 

V TT TT V 2 




3/2 r 


3 3 TT^ 

2-2G- - TT + (tt + 3) ln2 + - (ln2)2 - ^ ) = 0.0112528 , 
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with the Catalan constant 




(- 1 )’ 


-(2n + iy 

n=0 ^ ^ 


= 0.91596559 


The root approximant, respecting all given small-ce, as well as large-ce expansions, is 

X 5/6 \ 7/8 


I n /((I + + A^S'" + Ati^*A 




1/6 


with the parameters 


Co = 1.19055 , = 48.4532 , A 2 = 564.108 , 

/I 3 = 1088.39 , 2 I 4 = 1221.08 , A 5 = 796.791 , Ag = 256 . 

We estimate the accuracy of the root approximant comparing it with the numerical data from Ref. [38] 
and find that its maximal error is only 0.9%. Note that Fade approximants cannot be used in the case of 
harmonium, since the small-variable and large-variable asymptotic expansions are incompatible. 


6.2. Energy of two-electron spherium 

The two-electron spherium is a system consisting of two electrons that are confined to the surface of a 
sphere of radius R. The ground-state energy of the system [23,39] possesses the small-radius expansion 

E[R) ~ cq -f ciR -f C2R^ + c^R? {R —^ 0) , 

R 

in which 

Co = 41n2 - 3 = -0.22741128 , Ci = 8(ln2)^ - 401n2 + 24 = 0.11773689 , 

C 2 = -0.05027560 , C 3 = 0.01395783 . 

The coefficients C 2 and C 3 can also be expressed in closed forms that, however, are too much cumbersome 
[39], because of which we give here only their numerical values. 

In the large-radius limit, the energy has the expansion 

2 WR “ ^ “ 128i?5/2 (R^oo) . 

The root approximant can be writen in the form 

~ ^ + *70 ^ + A^R^^ -f A^R^ 

where 

= 1.05188915 , A 2 = 0.56453530 , Tig = 0.36000617 , 

TI 4 = 0.12606787 , TI 5 = 0.01946301 . 

Comparing this expression with numerical data [39], we find that the maximal error occurs at i? = 20, 
being only 0.1%. The best Fade approximant P^/^{-\fR) is much less accurate, having the maximal error, 
also at i? = 20, but an order larger, 1.5%. 
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7. Discussion 

We have deseribed a simple and general method for interpolating funetions between their 
small-variable and large-variable asymptotie expansions. The method is based on the eonstruetion of 
self-similar root approximants enjoying the general form 

fki^) = fo{x) (((... (1 -f AixY^ + A2X^Y^ + A'iX^Y^ -f ... -f AkJ^y^^ . 

All parameters Ai ean be uniquely defined through the eorresponding asymptotie expansions. By 
ehanging the variable, it is easy to inverse the expansions between the small-variable and large-variable 
limits. 

Our aim has been to suggest a method that would involve no fitting parameters. This is espeeially 
important in those eomplieated eases, where numerieal data in the whole region of the variable are 
not available. The absenee of fitting parameters makes our aproaeh different from other intrepolation 
methods, sueh as the Cioslowski method [23]. 

We have demonstrated the method of root approximants by several examples, whose strueture is 
typieal for many applieations, ineluding the hard-eore seattering problem, Debye funetion, Fermi-Dirae 
integral, Fekete-Szego problem, zero-dimensional oseillator, homogeneous eleetron gas, harmonium 
atom, and spherium. 

We have analyzed several more problems, e.g., the interpolation of the polaron mass between 
weak-eoupling and strong-eoupling limits studied earlier by the Feynman variational proeedure [40] 
and by other methods [41-44]. Our approaeh provides approximations, whose aeeuraey is eomparable 
or better than that of other methods, being at the same time more simple. 

Generally, the suggested method provides the aeeuraey not worse than the method of Fade 
approximants and in the majority of oases is more aoourate than the latter. 

Exoept the root approximants of the general form (9), we also have oonsidered additive approximants 
represented by the sums 

(M+N)/2 

fM/N^.^) ~ Aiil + . 

i=l 

This type of expressions ean be oonsidered either as additive root approximants or an additive variant 
resulting from self-similar faetor approximants [45] . 

For example, in the ease of one-dimensional eleetron gas, the oorrelation energy is approximated as 

^2/2(g) = ^i(l +^-f ^ 2 ( 1 -f , 


with the parameters 

= -0.044941 , A 2 = 0.017526 , = 0.266023 , B 2 = 0.133344 . 

This expression has the maximal error of 11%. However a more detailed analysis of sueh additive 
approximants requires a separate investigation, whieh is out of the seope of the present paper. 
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